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Abstract. The possibility to simulate lepton number violating supersymmetric models has been introduced 
into the recently updated Pythia event generator, now containing 1278 decay channels of SUSY particles 
into SM particles via lepton number violating interactions. This generator has been used in combination 
with the AtlFast detector simulation to study the impact of lepton number violation ) on event 
topologies in the ATLAS detector, and trigger menus designed for _^-SUSY are proposed based on very 
general considerations. In addition, a rather preliminary analysis is presented on the possibility for ATLAS 
to observe a signal above the background in several mSUGRA scenarios, using a combination of primitive 
cuts and neural networks to optimize the discriminating power between signal and background events over 
regions of parameter space rather than at individual points. It is found that a 5cr discovery is possible 
roughly for mi/2 < 1 TeV and mo < 2TeV with an integrated luminosity of 30fb~^, corresponding to one 
year of data taking with the LHC running at "mid-luminosity", £ = 3 x 10''^ cm~^s~^. 
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1 Introduction 



Among the primary physics goals of the Large Hadron 
Collider at CERN, scheduled to turn on some time during 
2006, are the exploration of the SM Higgs and top quark 
sectors, and the direct exploration of the TeV scale with 
the emphasis on Supersymmetry. 

The motivations for believing that Supersymmetry is 
indeed a property of Nature are many, most importantly 
the natural and exhaustive extension of the Poincare group 
furnished by the Supersymmetric operators the possi- 
bility for exact unification of the SM gauge couplings as 
required by GUTs, the requirement of Supersymmetry for 
anomaly-cancellation in string theories, and perhaps most 
importantly, its providing a natural and elegant solution 
to the problem of scale hierarchy in the SM. 

For this latter problem to be solved without unnatu- 
ral finetuning order by order in perturbation theory, the 
sparticle masses must lie at or below the TeV scale, and 
so there is ample reason to believe that SUSY should be 
observable at the LHC (see e.g. Q). 

The potential for hadron colliders to observe Super- 

fmmetry has been studied in great detail (see e.g. 
), yet studies tend to concentrate on the simple MSSM 
framework where both baryon and lepton number conser- 
vation are imposed. In section]^ we summarize the theoret- 
ical situation, pointing out that it is by no means obvious 
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that baryon and lepton number should be conserved in 
supersymmetric theories. 

The possibility to study lepton and baryon number 
violating Supersymmetry has already been included ^ in 
the Herwig event generator 1^ and in SusyGen 0, and 
several studies exist for run II of the Tevatron (for reviews, 
see 1^,^, latest experimental result s px|') and to a lesser 
extent for the LHC (see e.g. ]ll[0|l||j). 

However, it has not before been possible to study the 
full range of decays in these supersymmetric models us- 
ing the Pythia framework (and thus the string fragmen- 
tation model), and no dedicated LHC studies have been 
performed regarding triggers and sensitivity for these the- 
ories beyond cases where the LSP was simply forced to 
decay. In this paper, we go to the general case where all 
sparticles (excepting the gluino) are allowed to decay via 
lepton number violating couplings, introducing 1278 lep- 
ton number violating decay channels of sparticles to par- 
ticles in the Pythia event generator. 

For the purpose of studying the experimental signa- 
tures arising from the lepton number violating decay chan- 
nels, and to quantify to some extent the "discovery poten- 
tial" of the LHC, some selection of benchmark points is 
needed. 5 mSUGRA points and, for each of these, 9 sce- 
narios for the 36 coupling strengths are presented and 
discussed in section ^. 

From this point on, the analysis relies on technical as- 
pects of the detector design (or rather, their representation 
in crude simulation algorithms), and so we are forced to 
distinguish between the ATLAS and CMS experiments. 
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Assuming to a first approximation that these experiments 
will have similar capabilities, we choose to concentrate on 
ATLAS in the remainder of the paper. 

A first step is to define the data sample inside which 
-^-SUSY should be searched for. To this purpose, a selec- 
tion of triggers designed for ^I-SUSY scenarios are pro- 
posed for "mid-luminosity" running of the LHC (£ — 
3 X 10'^^ cm~^s~^) in section ||. 

Before going further, we note that the conventional 
procedure of studying a few benchmark points in detail 
is insufficient in the enlarged parameter space opened up 
by the L- violating couplings, even more so since the pres- 
ence of such couplings also remove important cosmological 
contraints on the mSUGRA parameters. In particular, the 
Lightest Supersymmctric Particle (LSP) is no longer re- 
quired to be neutral | p5[ , and essentially no bound on the 
mSUGRA parameters can be obtained by requiring a relic 
density less than or equal to the density of dark matter in 
the universe. 

In section ^ we make a first attempt at coming to 
grips with the size of the parameter space, proposing a 
method which relies on neural networks and the grouping 
of individual scenarios into classes, whereby regions rather 
than individual points become the studied objects in pa- 
rameter space. In this work, we define just 2 classes of - 
SUSY models, each containing 15 individual scenarios and 
a class consisting of 5 MSSM scenarios for reference. For 
each class, a neural network is trained with post-trigger 
events from all the scenarios in the class against back- 
ground events, allowing the network to pick out general 
qualities common to each class without over-fitting to a 
particular model. Although we study only a few scenar- 
ios here, the usefulness of this method should become ap- 
parent when considering the requirements posed by more 
comprehensive scans over parameter space. 

After some remarks pertaining to the dangers of us- 
ing neural networks, we present results for the estimated 
ATLAS sensitivity, S'/\/5 + B, for all scenarios with an 
assumed integrated luminosity of 30fb~^. The subsequent 
step of pinning down the model parameters once a signal 
has been observed, is clearly more model dependent and 
is not covered here. An up-to-date review can be found in 

i- . 

A brief outlook and concluding remarks are given in 
section H. 



2 L-Violating SUSY 

In this part of the paper, we do not present new results. 
Rather, we summarize some theoretical considerations con- 
cerning i?-violation which, although well documented in 
the literature, may not be in the active memory of all 
readers. 

Despite the many attractive features of Supersymme- 
try, most notably that it provides a natural solution to the 
hierarchy problem, it is well known that the most general 
supersymmetric Lagrangian (containing all terms obeying 
SU{3)c X SU{2)l X U{1)y gauge invariance and Super- 
symmetry) is utterly incompatible with experiment, re- 



gardless of whether the assumption of minimal particle 
content is made or not. 

The reason is that baryon and lepton number conser- 
vation is not guaranteed by any of the symmetries just 
mentioned, and the accidental conservation of these quan- 
tum numbers in the SM does not hold when extending 
the SM with Supersymmetry; the full SUSY Lagrangian 
contains renormalizable lepton and baryon number violat- 
ing operators which are suppressed only by the SUSY 
breaking mass scale squared. Unless a much more powerful 
suppression mechanism is also at work, these operators re- 
sult in a lifetime for the proton which would be measured 
in fractions of a second, to be compared with the experi- 
mental bound Tp > 1.6 x lO'^^yr (at 90% CL) obtained by 
Super-Kamiokande in the p ~* e+7r*' channel p^ . 

The possibility that the couplings responsible for pro- 
ton decay are just naturally small but not zero is almost 
out of the question. Assuming Msusy ~ 1 TeV, the prod- 
uct of the ^ and _^ couplings involved in proton decay is 
required to be less than 10"^*^ |jl^. Since at least one of 
the couplings is then forced to have a value below 10~^^ 
without any obvious suppression mechanism at work, it 
is more natural to assume that there is some additional 
symmetry in the theory, giving zero couplings either for 
the i?-violating terms, the L-violating terms, or both. 

Thus, in the Minimal Supersymmetric Standal Model 
(MSSM) , a discrete symmetry, i?-parity ||l^ , is customar- 
ily imposed which ensures the conservation of baryon as 
well as lepton number in the supersymmetric Lagrangian. 
Since there are strong indications that non-gauged sym- 
metries are maximally violated by quantum gravity effects 
(leading to wormhole- induced proton decay |^) whereas 
gauged symmetries are totally stable against such effects 
[ pi| , we shall here assume that the proton-protecting sym- 
metry is in fact a so-called discrete gauge symmetry [ p2| . 
In this case, it is of importance to note that although im- 
parity conservation leads to a comparatively simple phe- 
nomenology and a natural dark matter candidate (the 
LSP), the problem with proton decay is not satisfacto- 
rily solved when embedding the supersymmetric theory 
into more fundamental frameworks (such as GUTs) where 
baryon and lepton number violation can appear at some 
higher scale. 

In an effective theory valid around the electroweak en- 
ergy scale, the interactions mediated by super-heavy res- 
onances associated with a higher-scale theory, take the 
form of non-renormalizable operators (i.e. operators of 
mass dimension d > 5) suppressed by d — 4 powers of 
the high scale. Such operators, violating baryon and/or 
lepton number are generally not forbidden by i?-parity 
conservation (see fig. and it has been demonstrated 
that operators of dimension 5, being suppressed by only 
one power of the high (e.g. GUT) scale, will cause too 
rapid proton decay unless their couplings are suppressed 
by several orders of magnitude . This should, by itself, 
provide a powerful argument for exploring alternatives to 
i?-parity. 

In Q, a systematic analysis of discrete gauge sym- 
metries was carried out with the surprising conclusion 
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Fig. 1. An example of proton decay proceeding via a 5- 
dimensional operator violating both L and B but conserving 
R. Due to the Majorana nature of the neutralino, it is drawn 
without an arrow. 



that only i?-parity and a symmetry equivalent to L- 
conservation are anomaly-free, although some model de- 
pendence of this analysis ||2^ implies that also conserva- 
tion of baryon number could be assured by viable sym- 
metries. We shall here denote the two possibilities by L- 
parity and i?-parity, though this should not be taken to 
mean that these symmetries stand for exact lepton number 
conservation and exact baryon number conservation (for 
obvious reasons to do with the matter-antimatter asym- 
metry of the universe). These names only reflect that in 
the supersymmetric part of the theory, the correspond- 
ing quantities are conserved. They may still be violated 
by GUT scale operators and will certainly be so by elec- 
troweak sphalerons. By far the most interesting observa- 
tion, however, is that both L-parity and i3-parity but not 
i?-parity can forbid the dangerous dimension 5 contribu- 
tions to proton decay. 

Since a large number of current GUT and Planck scale 
theories do contain baryon and lepton number violation 
exterior to the supersymmetric MSSM framework, which 
will result in exactly this kind of contributions, it is im- 
portant that experiments planning to explore SUSY are 
well prepared for the very distinct experimental signatures 
which are the hallmarks of i?-parity violating scenarios, 
most importantly the consequences of a decaying LSP, of 
single sparticle production, and of the removal of the cos- 
mological constraint that the LSP be neutral. 

In this article, we assume the framework of the MSSM 
with the modification that lepton number is violated. This 
leads to the MSSM interaction Lagrangian being enlarged 
by the following terms [|l8j : 

-C^ = Xijk {^Li^L]e*Rk + ^ueLjeRk + I'LieLjeRk - i ^ j) 

'^Kjk {yLidLjd*Rk + VLidhjdRk + VLidLjdRk 

-^Ri'^Ljd*Rk - SLiULjdRk - e-LiULjdRk^ 

+h.c. (1) 

where A: are generation indices (summation implied) 
and \ijk is antisymmetric in its first two indices. The 
terms in the first line of the above equation are custom- 
arily denoted LLE terms and the terms in second LQD, 
in reference to the superfields appearing in the superpo- 
tential from which eq. (|]) is derived. All of these terms 
contain two (SM) fermions and one (SUSY) scalar. The 
oddness under R is thus directly visible here, since the SM 



fields have i? = -f 1 and the SUSY fields R =-l. Also, it 
is clear that each sfermion can decay in a number of ways 
to two SM fermions via these couplings. L- violating neu- 
tralino, chargino, and gluino decays are forced to proceed 
via one or more intermediate scalar resonances. 

A complete list of decay modes and full matrix ele- 
ments for all decays^ of MSSM sparticles to SM particles 
can be found in p6| . Of these, we have implemented all 
but the gluino decays into the Pythia event generator, 
publically available from version 6.2, with documentation 
included in the recently revised Pythia manual |27 . The 



conclusions of the present paper are not expected to be 
changed significantly by the inclusion of IJj gluino decays 
since the gluino is typically heavy and thus has a number 
of other, unsuppressed decay channels available. Details 
about the Pythia implementation can be found in [BSj. 

Note that the sfermion decay matrix elements in [p6| 
have been checked analytically by the author, and that 
the matrix elements listed there are not directly appli- 
cable in Pythia. Both Pythia and Herwig (Isawig) 
follow 1^,^ for the chargino and neutralino mixing con- 
ventions, but Herwig uses the opposite convention 
for the sfermion mixing angles, yielding a relative trans- 
position between the two programs. In addition, exten- 
sive counter-checks were made between the two programs 
to make sure they agree|^. Some non-negligible differences 
exist: 

1. Isawig does not use running masses in the evaluation 
of higgsino-type couplings whereas Pythia does. This 
has been observed to lead to substantial differences 
between the widths calculated by the two programs 
(in rare cases factors of 2 or more). 

2. Qfs at Mz is used by ISAWiG in the calculation of the 
gluino decays. In the forthcoming update of Pythia 
where gluino decays are included, as will be evaluated 
at the mass of the gluino. 

3. Pythia and Herwig (Isasusy) do not use the same 
RGE's, so many parameters can differ quite substan- 
tially at the EW scale from the same GUT input, cre- 
ating an "artificial" dissimilarity between the two pro- 
grams. Ideally, one should compare models with iden- 
tical EW scale mass spectra. 



3 mSUGRA Models 

The models used in this study have not been chosen among 
the ones initially suggested in the ATLAS Physics TDR 
This is partly due to the exclusion of most of these 
points by LEP (essentially from bounds on the Higgs mass), 
and partly since it is interesting to enable a direct com- 
parison between the capabilities of the LHC and other, fu- 
ture experiments. The 5 mSUGRA points shown in table ^ 
have therefore been selected among 14 points which were 
defined by the CLIC physics study group js^. More re- 
cently - too late to be included in the studies reported here 

^ This resulted in a few bug-fixes also in the Herwig code, 
so Isawig has had a few bugs in the ME calculation up to and 
including v. 1104. 
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- a new set of standard benchmarks points were proposed 
1^ which have now been adopted by both the CLIC and 
other hnear coUider communities. See 1^ for mass spectra 
analogous to table [|. To aid comparison, we briefly list 
the main differences/similiarities between our points and 
the new benchmark points. 

— Point 'A' is a (150 GeV) hghter version of P2. The total 
SUSY (pair production) cross section is approximately 
4 times larger than for P2: a a ~ 4ctp2. 

— Points 'B', 'C, T', and particularly 'G' are similar to 
Pg. upg ~ 0.5(tb w 5crc ~ 2(7/ ~ S.Sctg- 

— Point 'D' is dissimilar to all our points. It is closest 
to P2, but has a 500 GeV lighter gluino, and 400 GeV 
lighter squarks which increases the LHC cross section 
by an order of magnitude. 

— Point 'E' is a lower-mass version of F2, point 'F' a 
higher-mass version, ap^ ~ IStii? ~ O.ObcTE- 

— Point 'H' is similar to Pg but slightly heavier, giving a 
total LHC cross section[| an ~ j'^Pt 

— Point 'J' is a high-tan /3 version of P2 with very close 
to identical total cross section. The mass hierarchies 
are also very similar. 

— Point 'K' is a large-mi/2 sister of the only large-tan/? 
model included in this work, P12. The correspondingly 
larger masses yield a cross section uk ~ 0.04crpj2. The 
mass spectrum is similar to that of P7. 

— We have not included points at tan/3 — 45, but point 
'M' can be seen as an even heavier version of point 
'K'« P12 with a total cross section of only 0.1 femto- 
barns, and point 'L' has a mass spectrum not greatly 
different from Pg but with heavier squarks and gluinos 
giving a cross section an order of magnitude lower than 
for Pg. 

The total cross sections compared here include all MSSM 
pair production. Single sparticle production will give ad- 
ditional contributions depending on the strengths of the 
i?-violating couplings. 

One should keep in mind that all these points are de- 
fined for the MSSM, and as such have neutralino LSP's, a 
property which it has already been mentioned is not nec- 
essary in ^-SUSY scenarios. Though we shall not do so 
here, it is certainly advisable to explore the experimental 
consequences of non- neutralino LSP's in more detail. 

Both our points and the points in Q assume a van- 
ishing trilinear coupling at the GUT scale, i.e. Aq = 0. 
In connection with this work, a small study of the direct 
consequences of that assumption upon the results pre- 
sented here was performed. Varying ^0 between and 
500 GeV for P2, Pg, and P12 gave only a weak variation 
(0(5%)) of the semi-inclusive L- violating branching ratios 
(e.g. BR{xi qq'^ j), and so apart from the consequences 
of the distortion of the mass spectrum caused by Aq ^ 0, 
the main signatures (number of leptons, number of jets, 
etc.) should be only mildly affected by changes to this 
parameter. 

^ Point 'H' gave a LSP when using the GUT input pa- 
rameters in ||] with ISASUGRA. We therefore chose a 30 GeV 
larger value for mo, retrieving Xi as the LSP. 
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Pt 


P9 


P12 


F2 


CSUSY [fb] 


130 


3.9 
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110 
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xSOfb"^ 


3900 


114 
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Table 1. Selected points of analysis in the mSUGRA param- 
eter space, mass spectra as obtained with ISASUSY, and total 
SUSY pair production cross section at the LHC. 



In addition to the mSUGRA parameters come the 9 A 
(LLE) and the 27 A' (LQD) couplings for which the scenar- 
ios listed in table ^ have been studied for each mSUGRA 
point. The reason we do not consider couplings larger than 
10~^ is partly due to the present experimental bounds 
and partly due to the fact that some implicit approxima- 
tions in this analysis would break down for larger cou- 
plings: 1) the ^ couplings are not included in the RGE 
evolution of the SUSY parameters and 2) single sparticle 
production is not simulated. Note that the large- coupling 
(a) scenarios only observe the limits on individual cou- 
plings given in not the limits on products of the 
couplings. Thus, these scenarios are useful for studying 
the consequences of having many large couplings, but one 
should keep in mind that they are unrealistic in that not 
all of the couplings can be simultaneously large. 

On the other hand, if the couplings get signifi- 
cantly smaller than lO"**, the LSP lifetime can become so 
large that it decays outside the detector, mimicking the R- 
conserving scenarios which have already been extensively 
studied. For example, for the mSUGRA point P2, setting 



p. Z. Skands: Searching for L- Violating Supersymmetry at the LHC 



5 





Ifj c 


/~\ T T T* T T "XT 


M O D 11/ L a 






h 


n 






L L E 




\ ■ , 

k 




10-^ 




\f 















L Q D 















y 

^ijk 


10"^ 


10-^ 








L L E + L Q D 


^ijk 




10-^ 


sj rhe- ihe- ihe^ 


y 

"^ijk 


10^2 







Table 2. Selected points of analysis in the A — A' parame- 
ter space. The last column corresponds to the "natural cou- 
pling" scenario proposed in m = ^ = 126 CeV ^^'^ 
niq. = ^{rriuj + rndj)- These models (coloumn n) will be re- 
ferred to as "uLLE", "nLQD", and "uLLE -I- nLQD" in the 
text below. 

all A couplings to 10~^ and all A' couplings to zero results 
in a decay length for the LSP of re = 40 m. In the inter- 
mediate range, one may see the LSP decay directly inside 
the fiducial volume of the detector (see e.g. yet we 

abstain from relying on such a spectacular signature here 
so as not to be overly optimistic in our results. 

4 Triggers for / -SUSY 

A reasonable aim for the total -^-SUSY dedicated trigger 
rate is about IHz. We here focus on rates after the events 
have been filtered through the trigger system, i.e. we make 
no distinction between trigger levels. This is a technical 
issue which requires more detailed knowledge of the de- 
tector performance at mid-luminosity than is currently 
parametrized in AtlFast. Specifically, no parametriza- 
tion of the effects of pile-up at mid-luminosity is included, 
and so we here adopt a best-guess approach, performing 
the simulation without pile-up, and then multiplying the 
resulting trigger rates by a factor of 5/3 to estimate the 
true rates at £ = 3 x lO'^'^ cm^^s^^. This factor is based 
on the scaling exhibited by the inclusive electron, elec- 
tron/photon, and ^T+2jets trigger rates presented in 
chp.ll] from low to high luminosity. 

To retain as much generality as possible in the trig- 
ger definitions, it is sensible to use the information con- 
tained in the ^ superpotential terms rather than a se- 
lection of decay modes to define the trigger menus. The 
A (LLE) couplings are purely leptonic, and thus single 
sparticle production via these terms is not relevant for 
the LHC. Thus we are searching for at least two hard 
leptons (which, however, may be tauons), in most cases 
accompanied by I^t from escaping neutrinos. For the A' 
(LQD) terms, we expect at least two leptonic objects ac- 
companied by hard jets when the dominant production 
mechanism is MSSM pair production. Note that pure jet 
signatures would also be possible through single slepton 
production, (71^2 ~^ ^ — > (73^4. Single squark production 



EVENT RATES AND SAMPLE SIZES 
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a [mb] 


Rate [Hz] 


N 

^ ' gen 


QCD 2^2 








J,,, = 1 - 10 GcV 


55 


1.6x10** 


2.5x10^ 


= 10 - 75 GcV 


12 


3.7x10^ 


2.2x10^ 


pj, = 75 - 150 GoV 


5.5 xlO"^ 


1.7x10* 


1.4x10'^ 


PT > 150 GoV 


2.9x10""* 


8.7x10^ 


l.lxlO'^ 


ti 


6.2x10"'^ 


1.9 


5.9x10® 


z/w 


1.2x10"^ 


3.6x10^ 


1.8x10* 


zz/zw/ww 


1.2x10"'^ 


0.36 


5.9x10® 



Table 4. Numbers of generated events for the trigger study. 
The rates listed are total rates before trigger for £ — 3 x 
10^^ cm~^s~^. QCD events with the pr of the hard interaction 
below 1 GeV were not simulated. 



via LQD requries a lepton in the initial state and is thus 
suppressed at the LHC. 

The triggers so far investigated are listed in table ||, 
where background rates and efficiency ranges for all mod- 
els investigated are shown. For comparison, results for 
MSSM scenarios are also given. The event generation was 
performed with an augmented version of Pythia 6.155 
[|5| and AtlFast v. 2. 53 fs^. With respect to AtlFast, 
an attempt was made at obtaining more believable muon 
and electron reconstruction efficiencies by including by 
hand a muon efficiency of 95% independent of muon mo- 
mentum and an electron reconstruction efficiency of 80% 
for electrons with pT > 50 GeV and 70% for electrons with 
Pt < 50 GeV. These estimates are based on the ATLAS 
Physics TDR [|. 

Background cross sections, pre-trigger event rates, and 
sample sizes are listed in table ||. Though low-p^ QCD 
events have essentially no chance to pass triggers and 
much less the subsequent analysis cuts, the very high pu- 
rities required for SUSY signal extraction do not immedi- 
ately admit these events to be discounted entirely. Rather, 
a substantial sample of such events was generated with the 
object of placing an upper bound on the number of low-p^ 
QCD events remaining after cuts using Poisson statistics. 
This will be discussed in section ^, however some cau- 
tiousness should be employed in interpreting the bounds 
obtained, since the Monte Carlo is here being stretched 
far into the tails of its pt distributions. 

We do not consider triple gauge boson production. 
The cross section (excluding Higgs-induced production) 
is around 100 fb ||3^, i.e. at the same level as interesting 
SUSY cross sections, yet invariant mass cuts on pairs of 
jets would presumably be able to reduce this background 
considerably, and so we do not believe that these processes 
are dangerous as background sources. 

More detailed remarks and a large selection of plots 
of trigger rates and efficiencies versus thresholds can be 
found in ||. 

From the efficiencies in table ||, one easily sees how 
much cleaner the signatures of the purely leptonic (LLE) 
coupling are compared to the signatures involving quarks 
(LQD) where higher thresholds, due to the hadronic envi- 
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- 20 


% 


45 - 


- 70 


% 


10 - 


- 40 


% 
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Hz 
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% 
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% 
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Hz 
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- 80 


% 
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- 80 


% 
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% 


3j50 + mu20I 
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Hz 


5 - 


- 15 


% 


45 - 


- 60 


% 


12 - 


- 40 


% 


3j50 + e30I 


0.1 


Hz 


5 - 


- 10 


% 


15 - 


- 55 


% 


10 - 


- 35 


% 


3j75 + mel25 


0.1 


Hz 


30 - 


- 65 


% 


30 - 


- 70 


% 


40 - 


- 90 


% 


Total Rate 


2.1 Hz 


60 - 


- 90 % 


90 — 


99.9 


% 


60 - 


- 96 % 



/Combined Efficiency 

Table 3. Estimated trigger rates for background processes and trigger efficiency ranges for the various MSSM points and ^- 
SUSY scenarios studied for £ = 3 x 10"^^cm~^s~^. The nomenclature follows the ATLAS standard, where e.g. "mu451" means an 
isolated muon with pt > 45 GeV, "me" stands for missing (transverse) energy, and "3j50" means 3 jets, each with pr > 50 GeV. 
The total rate is smaller than the sum of the individual rates since there is a certain overlap, and the combined efficiencies can 
be larger than the individual efficiencies, since there is not a total overlap between the triggers. 



ronment, mean smaller efficiencies. From this we conclude 
that it would be of interest to extend the LQD study, 
examining whether 2-jet triggers and/or 3-object triggers 
could enhance the efficiency. 

Lastly, though the trigger proposals given here are de- 
signed exphcitly with ]// -SUSY in mind, they show a cer- 
tain overlap with triggers proposed for more conventional 
physics. The di-muon and 3 jets -I- lepton triggers, for ex- 
ample, have also been proposed for various Higgs searches. 
The di-electron trigger as well as 3 jets -I- electron are pro- 
posed to catch tt decays. Finally, the conventional SUSY 
searches also make use of multi- lepton, jets + I^t, and 
multi-jet signatures It is therefore not unlikely that 
the triggers proposed here can be incorporated to some 
extent into the conventional trigger programme. 



5 Discovery of ^ -SUSY at the LHC 

The main purpose here is to deliver an impression of what 
kind of signal strengths can be achieved at the LHC with 
30fb^^ integrated luminosity. To maintain generality in 
view of the more than thousand L-violating channels pos- 
sible, we do not discuss invariant mass reconstruction or 
measurements of the SUSY parameters in general. We fo- 
cus entirely on the isolation of candidate events through 
inclusive and kinematic event-shape variables in the at- 
tempt to obtain a statistically significant signal as com- 
pared with the background expectation after cuts. We do 
this for several models simultaneously (using neural net- 
works in the final step), an approach which is complemen- 
tary to the conventional one, where specific scenarios are 
studied one by one. 

Note that no attempt is made, beyond a crude worst- 
case estimate, to include the effects of pile-up in this anal- 
ysis. 



5.1 Missing Transverse Energy 

The post-trigger I^t distribution for the SM and its com- 
position is shown in figure ||a. Note that there are so few 
double gauge boson events that they are hardly visible on 
the plot. The peaks at I^t = 75 GeV and at I^t = 175 are 
due to the me75 and mel75 triggers becoming active. 

In figure ||b the distributions for the most "low-mass" 
mSUGRA point, Pg, are shown for P9(MSSM), Pga(LLE), 
and P9q(LQD). The degradation of the I^t signature in 
the ]/i models is evident, though one observes by com- 
paring with fig. ^ that a cut at low I^t values is still 
possible. Note that the I^t trigger peaks mentioned above 
are absent for the LLE scenario, this simply because the 
LLE scenarios do not rely to so great an extent on the I^t 
triggers, c.L table ||. 

The full range of mSUGRA models are plotted in fig- 
ure^. To enable the models to be shown on the same scale 
and to be distinguished from each other, the histograms 
have been normalized to unit area and smoothed. These 
plots are not intended to give detailed information, but 
rather to illustrate the spread between the models. Models 
with heavy squarks and gluinos (P7, F2, and to some ex- 
tent P12) have rather flat, MSSM-like signatures whereas 
the models with lighter sparticles are peaked at low I^t- 
This is an important point, since the light-sparticle sce- 
narios have high production cross sections and the heavy- 
sparticle ones low cross sections. A cut on I^t > 100 GeV 
seems a reasonable compromise between losing events in 
the peaked distributions (where we have many events any- 
way) and efficient background rejection required for the 
heavier scenarios (where we lose little by the cut but have 
fewer events). 



p. Z. Skands: Searching for L- Violating Supersymmetry at the LHC 



7 





Missing 

c) 

Fig. 2. a) and b): Ifr signatures for the SM and mSUGRA point 9a models (i.e. all relevant I/i couplings set to 10~^) normalized 
to 30fb-^ of data taking. "QCD LOW pr" means events from the 100 GeV < Pt < 150 GeV sample and "QCD HIGH pr" 
events from the pr > 150 GeV sample, c): Event distributions normalized to unit area for LLE, LQD, and the MSSM for all 
mSUGRA and _^ coupling points studied. The last row of plots is not intended to give detailed information, merely to illustrate 
the spread between the models. 



5.2 Leptons and Jets 

Due to the fact that the LSP decays, one expects an in- 
creased number of leptons and/or jets in the -SUSY sce- 
narios. In figures ||and ^the number of leptons is plotted 
versus the number of jets for each of the SM backgrounds 
and for P9(MSSM), Pga{LLE), and PgaCLQD). All events 
satisfy > 100 GeV. 

The absolute normalizations are, of course, very differ- 
ent for each of the SM plots. 

One notices that the number of reconstructed jets goes 
up to 15 on the plots in figs. || and IJ. The jet numbers are 
obtained using the AtlFast standard cone algorithm jet 
finder, employing a cone size of AR = y/ (Acf))'^ + {At])'^ = 
0.4 and with clusters up to pseudorapidities |?7| < 5 in- 
cluded only if their transverse energies are above EJp™ = 
10 GeV. Whether a resolution of up to 15 jets or more is 
possible (and reliable) in the finished experiment is ques- 
tionable, yet the important point remains that there is 
more hadronic activity associated with the LQD scenarios 
than is the case for the MSSM, and we should be able 
to distinguish between the two possibilities using any jet 
finder. 



The increased number of leptons in the LLE scenario as 
well as the increased number of jets (and a small increase 
in the number of leptons) in the LQD scenario relative to 
the MSSM are evident. Since the Z/W and high-pT QCD 
backgrounds have the highest cross sections, a cut remov- 
ing events in the lower left corner was performed, requiring 
^jots + Ni > 8 or Ni > 3. This cut has a rejection factor 
of about 20 for the SM events and efficiencies above 75% 
for all ^ -SUSY scenarios, again with a tradeoff between 
signal loss at low-mass points and efficient background re- 
jection for high-mass points. 



5.3 Additional Variables 

In addition to the above mentioned requirements, cuts on 
several other kinematical variables were performed before 
the neural network analysis below. We shall here only 
briefly define and comment these variables (see for 
details): 

— The Pt of the hardest object (jet or lepton) in the 
event was required to be greater than 200 GeV. 
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Fig. 3. Lepton versus jet multiphcity (see text) for back- 
ground events surviving the cut on Ifr- The numbering of the 
bins is such that the events with jets are in the bin to the 
right of the number and events with leptons are in the bin 
above the number 0. 



MSSM 



D D D I 
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Di 
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Fig. 4. Lepton versus jet multiplicity (see text) in mSUGRA 
P<j for events surviving the cut on Ipr- The numbering of the 
bins is such that the events with jets are in the bin to the 
right of the number and events with leptons are in the bin 
above the number 0. 



We define the "pr-weighted 4-object energy correla- 
tion" by: 



P4C 



1 f Ea E-i 



12 \E3 E2 



E2 

El 



(PTl + PT2 + PT3 + PTa) 

(2) 



where ...,i?4 (jpti, ■■■jPta) are the energies (trans- 
verse momenta) of the four hardest objects in the event. 
For pair-produced LSP's, one would expect four hard 
objects with more or less equal energies more often 
than would happen for background events. The en- 
ergy fractions are weighted by the average pt to give 
the variable some extra sensitivity against tt events. A 
cut at P4C > 200 GeV was placed. 
Events with Thrust greater than 0.85 were rejected. 
The Thrust calculation naturally suffers from the loss 
of particles down the beam-pipe at high pseudorapidi- 
ties. Only particles with \r]\ < 5 are defined as being 
inside the active calorimetry range in AtlFast, and 
only these particles are included in the Thrust calcu- 
lation. 

Events with Oblateness greater than 0.4 were rejected^. 
Events with Circularity less than 0.1 were rejected^. 



5.4 Neural Network Cuts 

It is clear that the above cuts do not define a dedicated 
search strategy. They were applied to all scenarios with 
no distinction between MSSM, LLE, and LQD models. A 
more comprehensive analysis would of course have to focus 
on each of these possibilities separately, and also some 
differentiation between the mSUGRA parameters would 
be required. In particular sparticles lighter than the top 
and ones heavier than the top could with advantage be 
searched for using separate strategies. 

This highlights the fact that the mSUGRA parame- 
ter space, including now the couplings, is not small. 
It was judged an uneconomical use of resources at this 
point to optimize a purely physics-based analysis for each 
of the many possibilities offered by this large parameter 
space. Instead, three neural networks were trained with 
post-trigger events to separate MSSM, LLE, and LQD 
scenarios, respectively, from the background distributions. 
Each network was thus trained not with one scenario but 
with equal numbers of events from 15 different scenar- 
ios with varying mSUGRA parameters and varying ]//- 
coupling strengths, allowing the networks to pick out gen- 
eral qualities common to each class without over-fitting to 
a particular model. 

This procedure is certainly prone to the dangers which 
are always present when using neural networks on simu- 
lated data, in that the input parameters (simulated ex- 
perimental observables) do not come with warning labels 
about where in their domains the simulated results can 



^ Note that Oblateness is calculated entirely in the transverse 
plane for hadron colliders. 

Same comment as for Oblateness. 
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be trusted and where not. However, exercizing caution 
on this point, the method proposed here could be a use- 
ful alternative for exploring theories with large parameter 
spaces. 

As inputs were used all of the above mentioned vari- 
ables along with the px of the four hardest jets and the two 
hardest leptons. The networks employed were single-layer 
perceptrons using the gradient descent learning algorithm, 
with biased linear input and output neurons, and hidden 
neurons activated according to a biased sigmoid ('logistic') 
function. At the end of each learning period, the networks 
were trimmed using the OBD (Optimal Brain Damage) 
prescription p9| . After ten such periods, the networks gen- 
erally showed a negligible difference in performance on the 
learning sample as compared with the performance on an 
independent monitor sample, indicating that over-fitting 
to the training samples is not a problem. 

The training samples consisted of 3000 post-trigger 
SUSY events, selected equally among the appropriate sce- 
narios, and 4000 background events, weighted according 
to their post-trigger cross sections to represent approxi- 
mately 10^ events in the learning algorithm. An output of 
zero (one) was the target for background (SUSY) events. 

The finished LLE (LQD) network showed a clear pref- 
erence for low-Thrust, multi-lepton (multi-jet), low-Oblate- 
ness events with a mixture of low and high IpT- For slightly 
higher Thrust values, low lepton momenta were preferred 
in the LLE case. For the LQD network, the most notewor- 
thy additional feature was a strong relaxation of all other 
considerations for very high jet energies. Note that the 
background levels at low Thrust values suffer from some 
theoretical uncertainty due to the fact that parton showers 
have been used to generate the multi-jet final states, an 
approach which in the past has been known to give too few 
multi-jet events. Though it is not clear that present gen- 
erators suffer from the same problems, it is important to 
note that there is a non-negligible and potentially danger- 
ous uncertainty in any analysis relying on parton-shower 
simulations for large jet numbers and low Thrust values. 

Processing now the full event samples remaining af- 
ter the cuts described above through each of the three 
(MSSM, LLE, and LQD) networks, a final cut was made 
requring outputs larger than 0.9. Since the learning sam- 
ples were small compared to the full event samples (e.g. 
about 100.000 SUSY events were generated per scenario 
and only a few hundred used for training) no effort was 
made to exclude those events which had participated in 
the learning process from the analysis. This is, in princi- 
ple, a source of error, yet we permit ourselves to ignore it 
due to the smallness of the learning samples and since any 
problems related to over-fitting have been minimized by 
OBD. For the SM backgrounds, the generated event sam- 
ples had been depleted considerably by the initial cuts, 
and only the tt and double gauge boson samples had any 
events remaining at all after the network cuts (note that 
the learning samples used for training - for both back- 
ground and SUSY - were post-trigger events). 

The procedure used to construct upper bounds on the 
actual event numbers was, for N events remaining in a 



sample, to calculate the mean, /i, of the Poisson distri- 
bution which has exactly 5% chance to result in N or 
fewer events. ^ is then interpreted as a "95% CL" upper 
bound on the number of events which could have passed 
the cut. For the low-pT QCD events, this number was then 
subjected to the same rejection factors under cuts as the 
high-pT sample (later the rejections for the tt sample were 
used for both), and the double gauge boson rejection fac- 
tors were used for the upper bound on single gauge events. 

Typically, for 30fb~^, between 500 and 1000 signal 
events remain after cuts for P27 Pi2, and F2. Many more 
of course remain for Pg because of the larger cross section, 
but there is no hope for P7 with only 114 events expected 
in total in 30fb~^ of data. Yet one should bear in mind 
that single sparticle production which has not been in- 
cluded here could significantly increase the cross section 
if the ]fi couplings are not much smaller than the gauge 
couplings. For a hadron machine like the LHC, this effect 
would only be big for the LQD terms since single slepton 
resonances would then be possible. Of course, if l/i terms 
are present, single squark production would be possible as 
well. 



5.5 Results 

We define the statistical significance with which a discov- 
ery can be made by 



P = 



S 



VS + B 



(3) 



and conservatively interpret P > 5 rather than the con- 
ventional S/ ^/B > 5 to mean that a 5a discovery will be 
possible. Using the event numbers obtained in the analysis 
(the estimates on B being 95% CL upper bounds) results 
in the P values listed in table |5[ 

In reality, P depends on (unknown) systematic QCD 
uncertainties (parton distributions etc) and should be cor- 
rected for the effects of pile-up, and so we can only be 
confident that a 5cr discovery is possible if P is somewhat 
larger than 5. Therefore, aside from working with the def- 
inition, eq. (^), we attempt to define a more pessimistic 
quantity in a very crude, ad hoc manner which we shall 

C0/H -Pcow ■ 

The non-inclusion of pile-up results in too optimistic 
estimates of S/B. To include an estimate of the reduction 
of this ratio, we rewrite eq. (||) to: 



(4) 



S 



V^TbJs 



where we now include the effects of pile-up by multiply- 
ing B/S hy some factor. That twice as many background 
events per signal event could be passing the analysis if 
pile-up was included seems a reasonably pessimistic guess. 

Furthermore, assuming that the intrinsic uncertainty 
on both S and B coming from uncertainties on QCD pa- 
rameters will, to a first approximation, work in the same 
direction and with a comparable magnitude for both B 
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and S, we expect that the denominator in the above for- 
mula is not affected by this uncertainty, and so we include 
the QCD-related uncertainties by reducing the number of 
signal events in the numerator by a factor of 1.5, believing 
this to be an adequate worst-case estimate. This yields the 
following form for the "corrected discovery potential" : 

Of course, this quantity should not be taken too seriously. 
We list it in table |^ merely to show the effects of the stated 
factors on the discovery potential, i.e. a reduction of S/B 
by a factor of 2 combined with a reduction of both S and 
B by 2/3. 

With regard to the sensitivity of this analysis on sys- 
tematic uncertainties on the normalization of B, note that 
only for Pg do we have an S/B 1. For the P2, P12, and 
F2 scenarios S/ B lies between one half and unity. We do 
not estimate this to be a serious problem since our back- 
ground estimate is only a 95% CL upper limit (larger event 
samples would most likely bring the high-p^ QCD compo- 
nent down) and since the background normalization can 
presumably be determined to better than a factor 2 (as is 
used in Pcorr) using complementary regions of phase space 
where the SUSY population is low or vanishing. 

What one can see in table I is again that the LLE 
scenarios typically yield much purer event samples than 
the LQD ones, due to the hadronic environment at the 
LHC. It is also worth noticing that it is not easy, from 
this analysis alone, to discriminate between MSSM, LLE, 
and LQD scenarios. This is due to the fact that the three 
networks were trained independently, i.e. they were taught 
to reject SM events but not events from the other scenario 
types. 

A suggestion for how we might have gained further 
insight would be to construct one network with one out- 
put node for each class of models, requiring an output of 
(0,0,...) for background events, (1,0,0,...) for events of class 
1, (0,1,0,0,...) for events of class 2, and so on. One would 
then be able to construct a measure for the relative prob- 
ability of the events surviving cuts belonging to a given 
class of models, by comparing their n-dimensional distri- 
bution in the network output space with the expected 
background distribution. 

6 Outlook and Conclusion 
6.1 Outlook: 

Though some preliminary studies have been performed in 
the present work, the field is large and many important 
tasks remain. Recently, a significant theoretical effort has 
been dedicated to studying the case of baryon number 
violation |^,^,^. Both single sparticle production and 
baryon number violating SUSY will be included in the 
Pythia generator in the foreseeable future. 

This will hopefully spur more LHC activity, with com- 
prehensive studies of both lepton and baryon number vi- 
olating scenarios. Preliminary studies indicate a lessening 



of the reach of the LHC in baryon number violating sce- 
narios , and it would be of interest to explore this reach 
with the full range of production and decay mechanisms 
included. In addition, studies should certainly be carried 
out for mSUGRA benchmark points without neutralino 
LSP's. 

Finally, triggers dedicated for ^-SUSY should be in- 
corporated in the standard ATLAS trigger programme. 

6.2 Conclusion: 

1278 decay modes of Supersymmetric particles into Stan- 
dard Model particles through lepton number violating cou- 
plings have been implemented in the Pythia event gener- 
ator. Combining this augmented version of the generator 
with a crude simulation of the ATLAS detector, trigger 
menus for mid-luminosity running of the LHC have been 
proposed and seen to have a high acceptance of supersym- 
metric events in several L-violating SuperGravity scenar- 
ios while still giving event rates in the IHz region. 

Taking these trigger menus as basis, the possibility for 
a 5(T discovery after 30fb~^ data taking was estimated 
for each investigated model, using a technique based on 
simultaneous classification of groups of models rather than 
detailed studies of single models. 

For cross sections down to 10^^" mb (or, roughly, for 
points with mg < 2TeV and "011/2 < ITeV) it was found 
that a 5(7 discovery was possible for all scenarios with 
30fb~^ of data. It is estimated that uncertainties related 
to QCD parameters or pile-up in the detector, both of 
which have not been taken into account in the present 
analysis, could not significantly affect this conclusion. 

L-violating decays of the gluino and the possibility of 
invoking P-violation will be implemented in the Pythia 
generator within a foreseeable future. 
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10.5/6.5 


pLQD 

7n 


1.0/0.6 


1.0/0.6 


1.1/0.6 










pLQD 
-'^9a 


116/91.6 


153/122 


125/99.0 










pLQD 
-^96 


113/88.7 


151/121 


123/97.6 










pLQD 
^9n 


113/88.5 


131/104 


110/86.6 










pLQD 

-■120 


15.7/10.0 


19.2/12.5 


20.9/13.7 










pLQD 
^126 


15.8/10.1 


19.4/12.6 


21.1/13.9 










pLQD 
■'^12n 


16.6/10.6 


18.6/12.1 


20.9/13.8 










pLQD 
^2a 


7.0/4.2 


9.5/5.9 


10.5/6.5 










^2b 


7.0/4.2 


9.4/5.8 


10.5/6.5 










pLQD 


6.9/4.2 


8.8/5.4 


10.3/6.4 











Table 5. ATLAS discovery potential and corrected discovery potential (see text) for all SUSY scenarios investigated using each 
of the three networks. The numbers shown correspond to an integrated luminosity of 30fb~^. 
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